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Abstract 

Next-to-leading order corrections to fragmentation functions in a light-cone gauge 
are discussed. This gauge simplifies the calculation by eliminating many Feynman 
diagrams at the expense of introducing spurious poles in loop integrals. As an 
application, the short-distance coefficients for the color-octet 3 S± term in the frag- 
mentation function for a gluon to split into polarized heavy quarkonium states are 
ii ■ re-calculated to order a 2 s . We show that the ill-defined spurious poles cancel and 

the appropriate prescriptions for the remaining spurious poles can be determined 
by calculating a subset of the diagrams in the Feynman gauge. Our answer agrees 
with the recent calculation of Braaten and Lee in the Feynman gauge, but disagrees 
with another previous calculation. 
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1 Introduction 

Factorization theorems for inclusive single-hadron production [|IJ guarantee that the domi- 
nant mechanism for a hadron production with high pr is fragmentation , the production 
of a parton which subsequently decays into the hadron and other partons. Many of the 
theoretical uncertainties disappear in the high px region and fragmentation is, therefore, 
a nice probe of the hadron production mechanism. This process is described by a frag- 
mentation function D(z, fj,), where z is the longitudinal momentum fraction of the hadron 
and \x is a factorization scale. 

The earliest calculations of fragmentation functions exploited the fact that the frag- 
mentation function is independent of the production processes of the decaying parton. For 
example, the fragmentation functions for heavy quark[|J and for quarkonia @,§,0 were de- 
duced by comparing the production cross sections for the hadron with the form predicted 
by the factorization theorems for inclusive single-hadron production. A field-theoretical 
definition of the fragmentation function can be expressed as matrix elements of bilocal 
operators in a light-cone gauge J7| or, more generally, as matrix elements of nonlocal 
gauge-invariant operators Jl|. The Collins-Soper definition of the fragmentation function 
was first used by Ma to calculate a fragmentation function for heavy quarkonium in lead- 
ing order || and in next-to-leading order ||. The definition is particularly convenient 
for carrying out calculations beyond leading order in a s , and it allows for the calculation 
in the Feynman gauge. By using the Feynman gauge, one can avoid the problem caused 
by the ambiguity of the spurious pole of the gluon propagator in the light-cone gauge. 
On the other hand, one must calculate a number of diagrams which do not appear in the 
light-cone gauge. In higher-order corrections to a fragmentation function, the simplicity 
in the light-cone gauge is remarkable, provided the spurious poles are handled correctly. 
The spurious pole problem does not appear in tree-level real-parton corrections. 

However, when we calculate a virtual correction to an amplitude in the light-cone 
gauge, we should keep in mind the possible problem caused by the spurious pole. A naive 
way to evaluate such light-cone dependent integrals is known as the Cauchy principal 
value (PV) prescription. If we use the PV prescription, there exist the ill-defined spurious 
pole in some loop integral. Without introducing ad hoc assumptions, these integrals are 
not calculable using dimensional regularization. An elegant method has been presented by 
Mandelstam and Leibbrandt independently JlO|JTl||.0 The ML prescription made it possible 



to transform such an integral into a well-defined one. The derivation of the Altarelli-Parisi 
evolution of parton densities |15| is one of the best examples of the use of the light-cone 
gauge. All previous calculations beyond-leading order had to employ some prescription 
for the spurious pole. The PV prescription is used in Refs. |7| JT^Jl7| . The ML prescription 



is employed in the leading order |18|] and the next-to-leading order |19| , p0 . 

A safe way to fix the prescription for the spurious pole is to find the gauge transforma- 
tion relation from another gauge, such as the Feynman gauge, where all the poles are well 
defined. Fortunately, at least for the case of the fragmentation function calculation, one 
may make use of the gauge-invariant property of the fragmentation function of Collins 



Comprehensive reviews can be found in Refs. [|12|,|13lJl 



and Soper. One notices that the spurious pole in the light-cone gauge is transformed into 
the propagator of the eikonal operator in the Feynman gauge which makes the nonlocal 
operators of the decaying partons, quark or gluon, gauge- invariant. Therefore, it is not 
necessary to introduce a prescription for the spurious pole structure. Instead, by matching 
the result from the light-cone gauge, where the spurious pole has an ambiguity, with the 
well-defined Feynman gauge result, one can determine the structure by virtue of gauge 
invariance. 

In this work, we develop a way to calculate next-to-leading order corrections to a frag- 
mentation function in the light-cone gauge. Instead of employing any known prescription 
for the spurious pole, we determine the sign of ie in the spurious pole by using only gauge 
invariance. It is fixed by comparing the light-cone gauge result with that from another 
gauge where there is no such ambiguity. This method, which has not been used before, 
is derived straight-forwardly from the gauge independent definition of the fragmentation 
function given by Collins and Soper ]I|. By choosing the gauge-fixing vector n with van- 
ishing transverse component with respect to the momentum of the produced hadron, we 
observe that the ill-defined spurious poles in one-loop integrals disappear, at least in our 
example. The reason is that the transverse momenta of the final states are integrated 
out. As an application, we re-derive the short- distance coefficients for the color-octet 3 Si 
term in the fragmentation function for a gluon to split into polarized heavy quarkonium 
states to order a 2 . There are two previous calculations of this function which disagree 
with each other |9|,plf. We use both the Feynman gauge and the light-cone gauge. We fix 



the sign of ie in the spurious pole by matching the two calculations. Our results from the 
two gauges agree with each other before the evaluation of the loop integrals. We remove 
ultraviolet divergences using the MS renormalization procedure. Our result agrees with 
the recent calculation of Braaten and Lee [ET] . 

This paper is organized as follows. In section [|, we give a short description of the 
spurious pole problem in the light-cone gauge, evaluate the one-loop renormalization 
constants in the light-cone gauge and compare our results with those based on the PV 
and the ML prescriptions. In section [|, we describe the method based on the Collins- 
Soper definition for the calculation of a fragmentation function in the light-cone gauge. 
As an application of our method, we give in section |] the result for the color-octet 3 Si 
term in the fragmentation function for a gluon to split into heavy quarkonia. A discussion 
is presented in section || 

2 One-loop correction in the light-cone gauge 

The light-cone gauge is a physical gauge where the gluon field A^ 1 has vanishing light-cone 
projection 

A-n = 0, (1) 

where n is an arbitrary light-like vector (n 2 = 0) appearing in the gauge-fixing term in 
the QCD Lagrangian. In the light-cone gauge, in which the fragmentation function was 



originally defined J7|, the eikonal line as well as the ghost decouples from the gluon, since 
the coupling, proportional to n M , is orthogonal to the gluon propagator. One draw-back 
of the light-cone gauge in higher-order calculations is the existence of the spurious pole 
1/k ■ n in the gluon propagator 
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where k is the momentum of the gluon. In order to evaluate a loop integral depending on 
the spurious pole, one need to introduce some assumptions. Various prescriptions for the 
spurious pole have been proposed. As a naive way, the PV prescription assumes 



k ■ n 



+ 



k ■ n + ie k ■ n 



ie 



An elegant method is known as the ML prescription HlO| , |TT 
spurious pole is re-expressed as 



(3) 
In this prescription, the 
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where n is a conjugate light-like vector satisfying n 2 = 0. Its spatial components are 
opposite in sign to those of n. The merit of the ML prescription is that it allows us for a 
proper Wick rotation to evaluate the integral in Euclidean space. In the ML prescription, 
the reduction of tensor integrals into scalar integrals, and the evaluation of scalar integrals, 
can be very involved due to the introduction of one more light-like vector n. The reason 
is that the light-cone vector has a non-vanishing transverse component compared to that 
of the propagator momentum. 

Let us classify the integrals involving spurious poles into two classes in view of the 
naive PV prescription which does not include any other assumption but (^). First, there 
is an integral with the ambiguous spurious pole. This integral can be regularized us- 
ing dimensional regularization imposing the PV prescription @. The integral can be 
evaluated by using the ML prescription as well. Once we choose n arbitrarily, some 
loop-integrals depending on this spurious pole are ill-defined even under dimensional reg- 
ularization within the PV prescription. For example, there is a scalar integral in the gluon 
self-energy diagram: 
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(5) 



where k is loop momentum. If we use the PV prescription without any further assumption, 
the integral @ is proportional to the ill-defined gamma function l/r(0) even though we 
use the regularized dimension D = 4 — 2e. Let us call this by the integral having the 
ill-defined spurious pole. In the ML prescription, the integral @ becomes a well-defined 
one. 



However, it is dangerous to impose an ad hoc prescription for these spurious poles. 
Once we employ such a prescription, we must check renormalizability and unitarity, the 
cancelation of infrared (IR) divergence, case by case. Fortunately, in the case of the 
fragmentation function calculation, there is a possibility to fix the spurious pole structure 
by using gauge invariance, so we do not have to depend on a specific prescription. In 
covariant gauges such as the Feynman gauge, there is no such problem. The spurious pole 
from the gluon propagator in the light-cone gauge is transformed into the propagator of 
the eikonal line in the Feynman gauge. And the ie sign for the eikonal line propagator is 
well defined. If we calculate a gauge-invariant quantity in the two gauges and compare 
the two results, we may fix the spurious pole ambiguity in the light-cone gauge. 

As an example, let us consider the next-to-leading order calculation of a fragmentation 
function. The natural choice of the light-cone vector n for defining a fragmentation 
function is n = (1, 0j_, — 1), which has vanishing transverse components relative to the 
daughter-particle momentum p = {po,0±,Pn)- in general, any fragmentation function is 
expressed in terms of the scalar products among vectors p, n, and n = (1, 0j_, 1). Since the 
transverse dependence is integrated out in the fragmentation function and the daughter 
particle is on its mass shell, the conjugate light-cone vector n is no longer independent of p 
and n. Therefore, the fragmentation function in the light-cone gauge becomes dependent 
only on p and n even if we use the ML prescription. In this case, the quark wavefunction 
renormalization constant Zq = 1 + 5Zq and the gluon propagator correction factor IT in 
one-loop level are given by 
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(8) 
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where n/ is the number of light quark flavors. The divergences are regularized using 
dimensional regularization with spatial dimensions N = 3 — 2e. Note that p denotes 
the momentum of the gluon for II and twice that of the quark momentum for Zq. The 
superscripts F and LC are used for Feynman gauge and light-cone gauge, respectively. The 
scalar integrals Iab- are given in Appendix A. The constants (0) and (0) for the Feynman 
gauge have no light-cone dependent integral while the light-cone-gauge counterparts (fJD 
and (|8|) have such integrals as Ibcd and Iabc- The integral Iacd does not appear in 
(!)"(!)■ m the light-cone gauge, Iacd appears in the vertex correction factor (|i~8|) . In 



the Feynman gauge, the light-cone dependent integrals appear only in the term (|2T| ) 
which involves the gluon coupling to the eikonal line. One must be careful about the ie 
prescription for the light-cone dependent denominator C defined in ( |A.4| ). In the Feynman 
gauge calculation, the sign of e is fixed since the propagators of both the gluon and the 
eikonal line are well defined. But we do not assume that the same sign is also valid 



for the quantities (|]) and (|8|) in the light-cone gauge. The appropriate ie prescription 
in the light-cone gauge can be fixed by comparing results for a gauge-invariant quantity 
calculated in the two gauges. 

The light-cone dependent scalar integrals Iabc, Iacd and Ibcd appearing in the light- 
cone gauge correction factors, possess important features. First, they do not have the 
ill-defined spurious pole, even though they have ambiguous spurious poles. The integral 
(I) with the ill-defined spurious pole cancels in the gluon self energy correction factor 
I1 LC in (^), so the integrals in (^ and (|8]) can be dimensionally regularized within the 
naive PV prescription (H) as well as the ML prescription. The reason why the ill-defined 
spurious pole does not appear in the light-cone gauge calculations in @ and (H) is that 
we have chosen the light-like vector n with vanishing transverse components relative to 
the momentum of the daughter particle. Furthermore, their values are independent of the 
sign of the ie in the definition ( A.4 ). Since the integral is invariant under the inversion 



/ — > —I and Ixyc\p^-p — —Ixyc, where X and Y are A, B or D defined in Appendix A, 
it is trivial to obtain the relations : 

Iabc = I abc* i Iacd = Iacd, Ibcd = Ibc*d- (10) 

where the definition of the integral Ixyc* is the same as that of the integral Ixyc except 
for the fact that C = (p — l)-n + ie is replaced by its complex conjugate C* = (p — l)-n — ie. 
At least in this case, the values of the scalar integrals agree with those evaluated by using 
the PV prescription : 

Ixyc - > ~ (Ixyc + Ixyc*) ■ (11) 

The independence on the sign in front of the ie in C might be accidental. If we use 
the ML prescription, each light-cone-dependent integral in ([H]) has ultraviolet (UV) and 
IR structures which are different from those shown in Appendix A. Effectively, the ML 
prescription transforms a double pole into an IR pole and makes the integral satisfy 
naive power counting rules. Note that Zq and IT are gauge dependent. The values in 
the Feynman gauge agree with well-known ones that can be found, for example, in Ref. 
pl| . The result using the ML prescription is known only for the UV poles. We have 
full agreement in the UV poles if we use the ML prescription : the gluon propagator 
correction term is proportional to the QCD beta function as IT = (33 — 2n/)a s /(127reuv) 
TT|.|2"2"|| and SZq C = a s /(37reuv) P3fl- All of them are listed , for example, in Ref. WA . 



Since we will use gauge invariance to determine the spurious pole structure, we do not 
proceed with the prescription dependence further. A thorough study of the application 
of the ML prescription to this problem will be presented elsewhere |25| . 



3 Collins-Soper definition and the light-cone gauge 

The fragmentation function D g ^ H (z, n) gives the probability that a gluon produced in 
a hard-scattering process involving momentum transfer of order /i decays into a hadron 

6 



H carrying a fraction z of the gluon's longitudinal momentum. This function can be 
defined in terms of the matrix element of a bilocal operator involving two gluon field 
strengths in a light-cone gauge [0]. In Ref. [0, Collins and Soper introduced a gauge- 
invariant definition of the gluon fragmentation function that involves the matrix element 
of a nonlocal operator consisting of two gluon field strengths and eikonal operators. One 
advantage of this definition is that it avoids subtleties associated with products of singular 
distributions. The gauge-invariant definition is also advantageous for explicit perturbative 
calculations, because it allows the calculation of radiative corrections to be simplified by 
using the Feynman gauge. 

The gauge-invariant definition of Collins and Soper for the gluon fragmentation func- 
tion for splitting into a hadron H is 

D „(z u) = (~ 9 ^ Z f + °° dx~e- ik+x ~ 

x(0\G+^0)£\0-) cb V H(zk+fi±) £(x-) ba G+»(0 + ,x-,0 ± )\0) . (12) 

The operator 8(x~) in (O) is an eikonal operator that involves a path-ordered exponential 
of gluon field operators along a light-like path: 



£(x ) ba = Pexp 
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(13) 



where A^{x) is the matrix-valued gluon field in the adjoint representation: [A M (x)] ac = 
if A b (x). The operator Vh( p +, p± ) in (|12"D is a projection onto states that, in the asymp- 
totic future, contain a hadron H with momentum p = (p + ,p~ = {m 2 H + p 2 j_)/p + ,p±), 
where tuh is the mass of the hadron. The hard-scattering scale /i in ( |P2| ) can be identified 
with the renormalization scale of the nonlocal operator. The prefactor in the definition 
(O) has, therefore, been expressed as a function of the number of spatial dimensions 
N = 3 — 2e. This definition is particularly useful when we use dimensional regularization 
to regularize ultraviolet divergences. If the production process of the hadron H can be 
described by perturbation theory, one can use the definition ([12]) to calculate the fragmen- 
tation function D g _> H (z, p) as a power series in a s . In Ref. [ffj, complete sets of Feynman 
rules for this perturbative expansion for quark and gluon fragmentation functions are 
given. By inserting the eikonal operator (|13|) , the operator consisting of two gluon fields 
with different locations becomes gauge invariant. At higher order in a s , there are nu- 
merous diagrams which have gluons coupled to the eikonal lines. In the light-cone gauge, 
the contribution of the eikonal operator disappears since the gluon decouples from the 
eikonal line. Therefore, there is a great reduction in the number of Feynman diagrams. 
On the other hand, the spurious pole contribution of the gluon propagator appears in the 
light-cone gauge. However, the gauge invariance of this definition fll2|) provides the gauge 
transformation of the eikonal line contribution in the Feynman gauge into the spurious 
pole contribution in the light-cone gauge. By comparing the final results for the gauge- 
invariant quantity D 9 ^h(z, p) from the two gauges, the spurious pole coming from the 
gluon propagator in the light-cone gauge can be fixed unambiguously. 



4 Application 

One remarkable example of a fragmentation phenomenon is charmonium production at 
the Fermilab Tevatron. The production rate of a heavy quarkonium depends on the cross 
section of a heavy quark pair QQ with small relative momentum. In high-energy pp 
collisions, the gluon production rate is dominant and the inclusive production of a heavy 
quark pair QQ via subsequent decay of this almost on-shell gluon is enhanced by the gluon 
propagator |4||. Furthermore, at leading order in a s , such a QQ pair created by the virtual 
gluon is dominated by a color-octet 3 Si state PS| . The color-octet 3 Si contribution has 
particular phenomeno logical importance. Braaten and Yuan showed that, in the gluon 
fragmentation function for splitting into triplet P-wave states, the infrared divergence in 
the short- distance coefficient of the color-singlet matrix element (Oi( 3 Pj)) can be avoided 
by including the color-octet 3 Si term |2"7|] . The production rate of direct J/ip and ip' at 



large pt at the Tevatron [28] is explained by Braaten and Fleming by introducing this 
(C 8 ( 3 ^)) term || 



In the NRQCD factorization formalism ||29|| , the fragmentation function D(z,/j,) for a 



parton splitting a heavy quarkonium is expressed as a linear combination of NRQCD ma- 
trix elements, which can be regarded as phenomenological parameters. The corresponding 
short- distance factors depend on z and are calculable in perturbation theory. Most of the 
phenomenologically relevant short- distance factors have been calculated to leading order 
in a s . They all begin at order a 2 s or higherf], with the exception of the color-octet 3 S± 
term in the gluon fragmentation function, which begins at order a s . Since the color-octet 
3 Si term dominates the high-p^ gluon fragmentation phenomena in heavy quarkonium 
production, the next-to-leading order correction of order o? s to this term is particularly 
important. 

As an application, we consider the next-to-leading order correction to the color-octet 
3 Si gluon fragmentation function for heavy quarkonium H. Since there is a discrepancy 
between CDF data and the leading-order prediction of the prompt J/ip polarization at 
large px where the gluon fragmentation contribution is important pD|pi|J52]J53| , the full 
NLO calculation of the polarized heavy quarkonium production rate is needed too. Un- 
fortunately, there are two different results for the color-octet 3 S± term P,|2"T|. Therefore, 
it is worth while to calculate this important function in an independent way. Since both 
previous calculations employed the Feynman gauge, we shall present our results in the 
light-cone gauge. In order to determine the appropriate prescription for the spurious 
poles, we use the result from the Feynman gauge. By comparing the two intermediate 
results before the evaluation of the light-cone dependent integrals, we fix the sign of ie in 
the spurious pole in the light-cone gauge. 



We use the same conventions as those presented in Ref. ]pl]] . We do not reproduce 
the description on the theoretical background of the fragmentation function for heavy 
quarkonium production in NRQCD factorization formalism which is well explained in 



Ref. [21|. Based on the NRQCD factorization formalism [2£], the fragmentation function 



2 for the color-singlet 3 Si channel, the short-distance factor begins at order a 3 , 



is written in a factorized form [^TJ : 

D^ H {z) = [(N - l)d T (z) + d L (z)] (Of ( 8 5i)>, (14) 

where dx and di, are the short-distance coefficients for the transverse and longitudinal 
contributions and (Og( 3 Si)) is the color-octet 3 S*i matrix element defined in Ref. |2U . 



There is only one lowest-order diagram in both Feynman and light-cone gauge, which 
is shown in Fig. [TJ. The circles connected by the double pair of lines represent the nonlocal 
operator consisting of the gluon field strengths and the eikonal operators in the definition 
(|l2l). The momentum k = (k + ,k~ ,k±) flows into the circle on the left and out of the 
circle on the right. The cutting line represents the projection onto states which, in the 
asymptotic future, include a QQ pair with total momentum p = (zk + ,p 2 /(zk + ), 0j_). The 
appearance of the diagrams for both gauges is the same in this order, since the circle 
should emit a gluon. With the Feynman rules of Ref. |l| and following the method of 
extracting the short- distance coefficients of the fragmentation function in Ref. [2I|, we 



can read off the order-a s terms in the short-distance functions dr(z) and d^z) as 

,2c 
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We have neglected the relative momentum of the heavy quark in the QQ rest frame so 
that the invariant mass of the pair is p 2 = 4mg. The LO results (|I5]) and ([H]) agree with 



previous calculations in the Feynman gauge plpl ] 



The Feynman diagrams for the fragmentation function for g — > QQ at order a 2 consist 
of virtual corrections, for which the final state is QQ, and real-gluon corrections, for 
which the final state is QQg. The diagrams with virtual-gluon corrections to the left 
of the cutting line are shown in Fig. |2|. The black blob in Fig. 0(a) includes the vertex 
corrections and propagator corrections shown in Fig. [| In the Feynman gauge, only 
the diagram in Fig. [|(b) vanishes, because the gluon attached to the eikonal line gives a 
factor of n M . On the other hand, all the diagrams except for Fig. [|(a) vanish in the light- 



cone gauge. If we use the threshold-expansion method of Braaten and Chen [|35|, we can 
simplify the structure of the expression without employing the projection method. With 
the threshold expansion, we can keep the full structure of color and spin. Here we utilize 



the dimensionally regularized threshold expansion method of Braaten and Chen |3^j36| 



With the Dirac equation and the usual methods for reducing tensor integrals into scalar 
integrals, we factorize each virtual correction diagram into the leading order diagram in 
Fig. [I], times a multiplicative factor. In the light-cone gauge, the ghost decouples since its 
coupling to the gluon is orthogonal to the gluon propagator (0), so the gluon propagator 
correction factor shown in Fig. 3(d) does not have ghost contribution. 

The virtual corrections contribute only to the transverse short-distance function dr{z) 



defined in 21 : 
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where 6Zq and II are defined in (P)— (^) and A is the vertex correction factor. The 
contribution from the remaining diagrams shown in Fig. 2 (b)-(e), which have gluon 
couplings to the eikonal lines, is expressed as A. Their values are expressed in terms of 
one-loop scalar integrals: 
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The explicit value of the vertex correction factor A F in the Feynman gauge shown in (|19|) 
agrees with the result in Ref. 



21 



The UV dependence of the vertex correction factor 
A LC in the light-cone gauge shown in (|T8| ) agrees with the result using the ML prescription 
in Refs. [23,24| where only the UV contribution is given : Aq C = —5Zq C = —a s /(3ire\jy). 
The integral Iaad has a Coulomb singularity as well as a logarithmic IR divergence due to 
the exchange of a gluon between the on-shell heavy quark and anti-quark. Dimensional 
regularization puts power infrared divergence like the Coulomb singularity to zero, so 
only the logarithmic IR divergence remains in the integral Iaad- Then the integral is 
effectively expressed by I add via the equation I add = (N — 4) Iaad- It is important to 
notice that various correction factors in the Feynman and the light-cone gauge involve 
different combinations of the same scalar integrals. Straight-forward sums for both gauges 
produce a common result 
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Thus the non-vanishing contributions from the gluon coupling to the eikonal line in the 
Feynman gauge, A F , is simply distributed to other correction factors in the light-cone 
gauge via additional gluon propagator terms. 

Since gauge invariance holds for both the virtual and the real-gluon corrections sepa- 
rately, the equality of the virtual corrections in the Feynman and the light-cone gauge is 
a consequence of gauge invariance. As we commented in the previous section, the light- 
cone dependent integrals in the Feynman gauge result have no ambiguities form spurious 
poles. On the other hand, we have not fixed the sign of the ie in the spurious pole of the 
integrals which are obtained in the light-cone gauge. Since we have found exact agree- 
ment between the two results in the two gauges, we may simply use the values obtained 



from the Feynman gauge calculation. Note that the integral Iabc disappears in (|22|) , so 
the only light-cone dependent integrals that survive are Iacd and Ibcd- The values of 
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these integrals are independent of the sign of the it in the definition of C in ( |A.4| ). The 



expansion of (|22|) in e reproduces the result of Braaten and Lee |2T 

d Cvirtual) fz)=d (LO) 
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(23) 



where /3 = (33 — 2ri/)/6. 

The Feynman diagrams for the real-gluon corrections to the fragmentation function 
for g — > QQ can also be calculated in both gauges. We draw the 5 left-half diagrams only, 
which must be multiplied by their complex conjugates to give a total of 25 diagrams. 
The real-gluon correction is a tree-level calculation. Therefore, there is no spurious pole 
problem. In the Feynman gauge, all 25 diagrams contribute, while only 9 diagrams in 
the light-cone gauge. In the latter gauge, diagrams [|(a) and |(b) vanish. The real-gluon 
correction contribution is also gauge invariant. Employing either gauge, we reproduce the 
real correction contribution given in Ref f2lj before the phase-space integral is performed: 



4 eal) 



(z) 



Tca s [i 



d 



(real) , 




dx- 



(l-z)/z 



X A 



(24) 
(25) 



where t = (1 — z)(zx + z — 1), x = 2q ■ p/p , q is the final-state gluon momentum, and 
p is the QQ momentum. The final results for the real-gluon correction contribution of 
Braaten and Lee are straight-forwardly reproduced: 
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(26) 
(27) 



The infrared divergence cancels after summing the real and virtual correction con- 
tributions shown in (|23| ) and (p6|) . Employing the MS scheme, a s and the operator are 
renormalized as in Ref. |H]]. After renormalization, the final answers for dx{z) and di{z) 
of Braaten and Lee 01 are reproduced: 



dr(z,fj) 
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where the coefficient A(/i) is 



^<H ln 4 + ir)- 

and P gg (y) is the gluon splitting function: 

1-z 



71 



PM 



6 



+ 81n2 + 61n 2 2 



Po 



(29) 



1-z) 



: + Z(1-Z) + ^5(1-Z) 

z o 



(30) 



The transverse term dr(z) in (|28| ) still disagrees with that of Ma 0. Our final answer for 

function 

d L (z,[i) 



the longitudinal fragmentation function is obtained by setting e — > in ([27]): 



5m? 



(31) 



The longitudinal term, 6^(2:) agrees with that of Braaten and Lee |2l| as well as that of 
Beneke and Rothsteinf3"T||. The dependence on the spectroscopic state of the produced 



quarkonium of this fragmentation function can be found in Ref. [21 



5 Discussion 

We have shown how the light-cone gauge can be used to evaluate the perturbatively cal- 
culable parts of a fragmentation function. As an application, we tested our method by 
evaluating the next-to-leading order correction to the color-octet 3 Si term in the gluon 
fragmentation function. We reproduced the recent result of Braaten and Lee |^TJ which 
disagrees with that of Ma |J. The light-cone gauge considerably simplifies the calculation 
procedure for both the real and the virtual corrections. At least at the one-loop level, 
the spurious pole problem can be resolved. This problem does not appear in the real 
corrections, because they come from tree-level diagrams, but it does appear in the vir- 
tual corrections. The gauge-invariant definition of the fragmentation function of Collins 
and Soper allows us to fix the ambiguities from spurious poles in the light-cone gauge 
by comparing with the result obtained in the Feynman gauge. We reduced the virtual 
correction in the color-octet 3 S± fragmentation function in the light-cone gauge to a lin- 
ear combination of scalar integrals. After naive cancelations among the scalar integrals, 
ignoring the ambiguity from spurious poles, the correction reduces to scalar integrals that 
are independent of the sign of e in the denominator k ■ n + ie. Thus the PV prescription 
gives the correct answer. To see if the ML prescription also gives correct answer requires 
explicit calculations of the scalar integrals including IR and finite terms. As a byproduct, 
the renormalization constants in the light-cone gauge were obtained at one-loop level. 
Their UV dependencies agree with the previous calculations within the ML prescription. 
They might be useful for other calculations, such as the next-to-leading order corrections 



to other fragmentation functions [38 . 
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Were it not for the problem of the spurious pole, one could reduce a large amount of 
the intermediate calculation by using the light-cone gauge. If we choose the gauge-fixing 
vector n with vanishing transverse components with respect to the momentum of the 
produced hadron, the one-loop integrals are remarkably simplified compared to the case 
using the ML prescription. We were able to determine the spurious pole structure without 
depending on a specific prescription by comparing with the result in the Feynman gauge. 
Of course, if we have to repeat the entire calculation using the Feynman gauge to fix the 
spurious poles, the light-cone gauge will not save any labor. However, there are only a 
small number of integrals which have the spurious pole problem, and this provides a way 
to save labor. First calculate the full contribution in the light-cone gauge in terms of scalar 
integrals, without specifying any prescription for the spurious poles. Then, calculate in 
the Feynman gauge only those diagrams where the eikonal line couples with one or more 
gluons. By comparing the two results, we can determine the appropriate prescription for 
the spurious poles in the integrals. One restriction of this method is that it can only be 
applied to gauge-invariant quantities. 

The study of high-p^ fragmentation phenomena has a significant potential to refine 
our understanding of hadron physics. It avoids the nontrivial resummations that compli- 
cate theoretical predictions for low-pr hadron processes. The signature of fragmentation 
dominance in high px charmonium production has been observed in Run I of the Teva- 
tron. In Run II of the Tevatron, as well as at LHC and at future colliders, there will 
be much better statistics of high-pr heavy meson events. Quantitative predictions for 
quarkonium production at high p? will require next-to-leading order calculations of all 
the phenomenologically relevant fragmentation functions. The light-cone gauge may be a 
powerful tool for carrying out these calculations. 
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A Integral Table 

In this appendix, we present the explicit values of the integrals encountered in evaluating 



the virtual-gluon corrections. Most of them are presented in Ref. [21]. But we reproduce 
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here for completeness. These integrals have the form 

d N+1 l 1 



<AB- 



(A.l) 



(2tt) n + 1 AB ■ ■ ■ ' 
where the denominator AB ■ ■ ■ can be a product of 1, 2, 3, or 4 of the following factors: 



A = I 2 + ie, 

B = (l-p) 2 + ie = l 2 -2l-p + Am 2 Q + ie, 

C — (p — I) • n + ie, 

D = (l- p/2) 2 -m 2 + ie = l 2 -l-p + ie. 



(A.2) 
(A.3) 
(A.4) 
(A.5) 



The momentum p is that of a QQ pair with zero relative momentum (p 2 = 4mg) and n 
is light-like (n 2 = 0). The integrals I a and Ib vanish in dimensional regularization. By 
symmetry under p — > I — p, we have I ad — Ibd- Some of the integrals can be reduced to 
ones with fewer denominators by using the identity A + B — 2D = Aitlq: 



im 2 Q I ABD = 2(7 AD - I AB ) 
The independent integrals that need to be evaluated are therefore 
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(A.8) 
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(A.12) 
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(A. 14) 



The subscripts on the poles in e indicate whether the divergences are of ultraviolet or 
infrared origin. 
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Figure 1: Leading order Feynman diagram for g — > QQ. 
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Figure 2: The Feynman diagrams of order a 2 s for g — > QQ with QQ final states. There 
are additional contributions from the complex-conjugate diagrams. 
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Figure 3: One loop correction diagrams for g* — > QQ. 
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Figure 4: The Feynman diagrams of order a 2 s for g — > QQ with QQg final states. There 
are a total of 25 diagrams, but only the left halves of the diagrams are shown. 
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